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1.  INTRODUCTION 


In  May  1976,  Tropical  Storm  Olga  drifted  slowly  westward 
across  the  Philippines  and  essentially  stalled  some  75  n  mi  north- 
northwest  of  Subic  Bay,  one  of  the  major  U.S.  Navy  installations 
in  the  western  North  Pacific.  As  this  rather  weak  tropical  cyclone 
(maximum  winds  of  40  kt)  contributed  southwesterly  winds  into  Subic 
Bay,  the  vulnerability  of  vessels  within  the  harbor  became  evident. 
With  maximum  surface  winds  of  only  35-45  kt  observed  in  the  harbor, 
five  Navy  ships  sustained  light  to  moderate  hull  damage  and  became 
part  of  the  total  $1.4  million  damage  to  ship  and  harbor  facilities 
caused  by  TS  01 ga  . 

With  the  multitude  of  anchorages  and  piers ide  berths  avail¬ 
able,  a  major  problem  for  ship  commanders  and  other  decision  makers 
in  Subic  Bay  was  where  to  place  vessels  within  the  harbor  complex 
so  as  to  minimize  stress  on  ship  and  port  facilities.  Knowledge 
of  the  stream! ine-isotach  conditions  within  the  harbor  for  winds 
coming  from  the  vulnerable  directions  of  south  through  west  would 
certainly  have  aided  t r. e  uecision- making  process. 

The  mesoscale  nature  of  the  problem  is  further  complicated  by 
the  complex  topography  of  the  Subic  Bay  area  (Figure  1).  The  bay 
is  surrounded  by  mountainous  terrain  to  the  west,  northeast,  and 
southeast,  with  passes  through  the  mountains  to  the  ea s t- no r t hea s t 
and  northwest . 

In  order  to  simulate  in  a  laboratory  setting  the  surface 
streaml i ne-i sotach  distributions  throughout  the  harbor,  a  scale 
model  of  the  Subic  Bay  area  was  constructed  at  the  Fluid  Dynamics 
and  Diffusion  Laboratory  (FDDL)  at  Colorado  State  University. 

The  information  derived  from  this  experiment  was  compared  with 
synoptic  analyses  of  strong  southwesterly  wind  situations  observed 
at  Subic  Bay.  The  experiment  results  also  were  compared  indepen¬ 
dently  with  the  results  of  a  high  resolution  primitive  equation 
model  designed  for  the  Subic  Bay  area. 

This  report  presents  the  finding  of  these  independent 
analyses,  in  order  to  provide  guidance  to  decision  makers  who  must 
plan  for  tropical  cyclone  stress  conditions  in  the  Subic  Bay  area. 
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Figure  1.  Subic  Bay  area  on  the  west  coast  of  the 
Philippine  Islands  at  14°  15'  N,  120°  14'  E. 
Anchorages  generally  are  available  within  the  bay 
area  for  depths  greater  than  60  ft.  Terrain  topo¬ 
graphical  features  and  bottom  contours  shown  in 
feet.  Arrow  points  to  Cubi  Point  Runway. 


Figure  2.  Topographical  model  in  wind  tunnel,  viewed 
from  the  south.  Grid  squares  in  bay  area  are  four 
inches  on  a  side  and  relate  to  0.82  n  mi  squares. 


2.  LABORATORY  MODEL  AND  RESULTS 


2.1  LABORATORY  MODEL 

The  1:15,000  scale  model  of  the  Subic  Bay  area  used  in  the 
experiment  is  shown  in  Figure  2.  The  terrain  was  constructed 
from  styrofoam  sheets  cut  to  match  contour  lines  of  the  topography; 
surfaces  were  smoothed  with  clay  and  painted.  The  model  was 
studied  in  the  Fluid  Dynamics  and  Difusion  Laboratory's  environ¬ 
mental  wind  tunnel,  shown  in  Figure  3. 

The  study  included  flow  visualization  using  oil  films  and 
smoke  (titanium  tetrachloride),  as  well  as  hot-wire  anemometer 
measurements.  Rotation  of  the  scale  model  facilitated  development 
of  surface  flow  patterns,  mean  wind  speed  isotachs,  and  gustiness 
values  for  the  entire  Subic  Bay  area  for  incoming  wind  directions 
from  south  through  west  (15°  intervals).  For  more  details  on  the 
experimental  technique,  the  reader  is  referred  to  Woo  et  al  .  , 

1  978  . 

Values  of  mean  wind  speed  and  turbulence  intensity^  (related 
to  gustiness)  were  examined  at  four  scale  elevations  --  .06,  .12, 
.16  and  .22  inches,  representing  heights  in  the  real  atmosphere 
of  70,  150,  195  and  275  ft,  respectively  --  for  the  60  locations 
shown  in  Figure  4  . 

The  mean  wind  speed  and  maximum  gust  profiles  taken  upstream 
of  the  bay  entrance,  are  characteristic  of  the  boundary  layer 
upstream  of  the  model.  The  boundary  layer  thickness  5  was  2.5 
inches,  corresponding  to  a  prototype  of  3125  ft  in  the  real 
atmosphere.  The  relationship  to  examine  the  vertical  variation 
of  the  horizontal  wind  speed  can  be  written  as 


where  U  is  the  mean  wind  speed  at  level  z,  is  the  mean  wind 
speed  at  the  top  of  the  boundary  layer  at  height  6,  and  n  is  a 


^Turbulence  intensity,  UT,  is  related  to  mean  wind  speed,  U  and 
maximum  gust,  Ug,  by:  =  U(2UT  +  1). 
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Figure  4.  Measurement  locations  for  wind  speed  and 
turbulence  intensity.  Topographical  elevations 
shown  in  f t . 
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variable  delineating  the  variation  of  the  mean  wind  speed  in  the 
vertical.  As  indicated  in  Bates  (1978),  data  from  hurricanes  have 
shown  n  to  vary  from  0.14  to  0.30.  An  n  value  of  0.22,  which 
fits  well  within  the  range  of  observed  values,  was  found  for  the 
vertical  variation  of  wind  in  the  laboratory  experiment.  Figure  5 
presents  the  mean  wind  speed  and  maximum  gust  profiles  of  the 
approach  flow. 


2.2  RESULTS 
2.2.1  Mean  Wi nds 

Figure  6  shows  examples  of  the  streamline-isotach  distribu¬ 
tions  derived  for  incoming  wind  directions  of  180°,  225°,  and 
255°.  The  figures  in  Appendix  A  show  the  streamline-isotach 
distributions  for  the  seven  incoming  directions  examined.  The 
distributions  shown  in  Figure  6  are  for  the  lowest  level  of  isotach 

measurements  (i.e.,  simulation  of  70  ft),  with  an  approach  wind  of 
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30  kt.  Approach  wind  speeds  for  other  levels  can  be  obtained 
from  the  profiles  previously  provided  in  Figure  5. 

In  general,  the  wind  speed  is  uniform  inside  the  bay  for 
southerly  winds  (Figure  6a),  except  to  the  lee  of  land  areas.  The 
wind  characteristics  become  more  and  more  complicated  as  the  wind 
becomes  more  westerly,  because  of  the  mountains  and  valleys  to  the 
west  of  the  Subic  Bay  complex.  Winds  stronger  than  the  approach 
wind  are  also  observed  near  the  mouth  of  the  harbor  for  the  more 
southwesterly  wind  directions  (see  Figure  6b).  For  westerly  winds 
(Figure  6c),  the  flow  has  a  slight  northerly  component  in  the 
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The  approach  speed  at  70  ft  is  purely  a  function  of  the  vertical 
variation  of  wind  from  the  relationship  shown  in  Figure  5.  The 
value  at  70  ft  is  0.43  the  value  at  the  top  of  the  boundary  layer. 
For  simplicity,  this  is  set  at  30  kt  for  ease  in  the  comparisons 
with  the  synoptic  and  dynamic  model  results  to  be  discussed  later. 
This  value  could  be  adjusted  either  up  or  down  with  resulting  wind 
distributions  appropriately  adjusted.  For  example,  if  the  incoming 
wind  is  45  kt  rather  than  30  kt,  then  the  20  kt  isopleth  within 
the  harbor  would  become  a  30  kt  isopleth.  The  laboratory  experi¬ 
ment  was  not  performed  at  varied  incoming  wind  speeds,  but  only  at 
varied  directions. 
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Figure  5.  Mean  wind  speed  and  maximum  gust  profiles  of 
the  approach  flow.  U  is  the  mean  wind  speed  at  level  z, 

U  is  the  mean  wind  speed  at  the  top  of  the  boundary 
layer  at  height  5,  and  U  is  the  maximum  gust  wind  value. 
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Figure  6.  Streaml i ne- i sotach 
distributions  for  wind  tunnel 
laboratory  experiment  for  wind 
directions  from:  (a)  180°, 

(b)  225°,  and  (c)  255°.  Wind 
speed  values  in  kt  are  for  an 
elevation  of  70  ft;  incoming 
wind  speed  is  30  kt. 


255°/30  kt 


upper  section  of  the  harbor.  The  mountains  generally  provide 
good  sheltering  for  the  region  near  the  west  of  the  bay  when  the 
wind  approaches  from  behind  the  mountains,  but  this  effect  is  not 
observed  on  the  east  side  of  the  bay.  In  some  cases  winds  even 
higher  than  the  approach  wind  are  observed  due  to  converging 
secondary  flows  through  the  valleys  between  mountain  peaks. 

2.2.2  Maximum  Gustiness. 

As  noted  earlier,  another  application  of  the  model  results  was 
the  development  of  the  capability  to  derive  maximum  gustiness  values 
incorporating  the  mean  wind  and  turbulence  intensity  measurements. 
For  example.  Figure  7  depicts  (a)  the  mean  wind  of  an  approach 
flow  of  255°  at  30  kt  and  (b)  the  derived  maximum  gust  values 
(highest  wind  over  a  one-minute  period).  Maximum  gusts  over  50  kt 
were  observed  in  the  central  bay  region.  In  the  sheltered  western 
portion  of  the  harbor,  maximum  gust  values  were  approximately 
twice  the  observed  mean  wind  values.  Near  the  harbor  mouth,  which 
is  much  less  influenced  by  terrain  than  the  interior  of  the  Bay, 
the  gust  factor  (gust  value)  was  approximately  1.4  the  mean  wind 
va 1 ue  . 

These  figures  can  be  compared  with  a  tropical  cyclone  gust 
factor  of  1.1  to  1.8  found  by  Gentry  (  1  953)  for  the  ratio  of  gust 
speed  to  five-minute  mean  speed.  Shiotani  (1975)  found  values 
ranging  from  1.2  to  1.5  for  a  typhoon  passage  as  observed  by  a 
meteorological  tower.  Atkinson  and  Holliday  (1977)  deduced  a 
value  of  1.25  for  a  maximum  sustained  40  kt,  one-minute  wind 
speed  for  a  tropical  cyclone  over  water. 

Ship  commanders  will  sometimes  find  it  more  important  to 
have  knowledge  of  the  gust i ness  associated  with  wind  conditions 
than  simply  knowledge  of  a  consistent  strong  wind.  Commanders  can 
adjust  to  consistent  strong  winds,  but  the  variable  forces  of 
gustiness  can  be  hazardous  to  a  ship  at  anchor  or  at  pierside. 
Analyses  of  the  maximum  gusts  for  the  seven  wind  directions 
examined  in  the  laboratory  experiment  are  shown  in  Appendix  B. 
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Figure  7.  Streaml i ne- i sotach  distributions  for  (a)  an  incoming 
wind  from  255°  at  30  kt  and  (b)  the  maximum  gust  values  in  kt 
Both  distributions  are  for  the  70  ft  elevation. 
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Figure  8.  The  vertical 
variation  of  wind  at  a 
point  1/3  n  mi  to  the 
east-northeast  of  Cubi 
Point  Runway,  based  on 
winds  from  south  through 
west.  The  mean  and  extremes 
of  the  differences  of  the 
wind  speed  at  any  elevation 
Z,  minus  the  wind  at  70  ft 
are  plotted  for  the  vertical 
variation  from  70  to  275  ft. 


Wind  Speed -Wind  Speed  |kt) 
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2.2.3  Vertical  Variation 


The  vertical  variation  of  wind  and  the  turbulence  associated 
with  wind  are  important  considerations  in  the  structural  design 
of  large  buildings  and  bridges  in  coastal  regions  and  harbors 
likely  to  be  affected  by  tropical  cyclones.  Personnel  at  coastal 
airport  locations  also  must  be  aware  of  the  vertical  wind  shear 
in  the  critical  lowest  few  hundred  feet  under  strong  wind  or 
typnoon  stress  conditions. 

In  the  Subic  Bay  model  experiment,  the  vertical  variation  of 
the  wind  throughout  the  harbor  for  each  wind  direction  was  deduced 
for  the  simulated  boundary  layer  from  70  ft  to  275  ft.  An  example 
of  the  vertical  variation  of  the  wind  for  one  particular  location 
(Point  A  in  Figure  4)  is  shown  in  Figure  8.  Assuming  south¬ 
westerly  winds,  this  location  is  downwind  approximately  1/3  mi 
from  the  end  of  Cubi  Point  runway. 
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3.  SYNOPTIC  ANALYSIS 


In  the  wake  of  the  damaging  southwest  winds  generated  by 
Tropical  Storm  Olga  in  May  1976,  the  Naval  Weather  Service 
Environmental  Detachment  at  Cubi  Point  coordinated  a  data  collec¬ 
tion  program  within  the  Subic  Bay  area.  The  aim  was  to  obtain 
extensive  synoptic  information  from  ships  located  throughout  the 
harbor  complex,  whether  at  anchor  or  at  pierside,  during  strong 
southwest  wind  conditions.  During  the  period  May  1976  through 
September  1977,  five  synoptic  situations  resulted  in  strong  south¬ 
westerly  winds:  four  were  associated  with  tropical  cyclones  in 
close  proximity  (Tropical  Storms  Olga,  21-25  May  1976;  Iris 
14-16  September  1976;  Sarah,  12  July  1977;  and  Dinah,  16  September 
1977);  and  the  fifth  included  a  period  of  unusually  strong 
monsoonal  southwesterly  winds  (17  May  1977). 

The  data  set  included  observations  from  ships  as  well  as 
information  from  a  recording  weather  station  located  near  the 
airport  runway  at  a  height  of  approximately  70  ft  (see  Figure  9). 
Since  ship  anemometers  vary  considerably  in  elevation  (for  example, 
some  carriers  have  their  anemometers  over  150  ft  above  mean  sea 
level),  all  observations  were  normalized  to  the  70  ft  elevation  by 
the  vertical  variation  of  tropical  cyclones  discussed  by  Bates 
(  1  978  ).  Figure  1  o'  shows  the  normalized  profiles  for  the  vertical 
wind  variation  for  tropical  cyclones  over  land  and  over  water  as 
described  by  Bates.  A  mean  value  between  these  two  curves  was 
used  in  this  segment  of  the  study.  The  normalized  profile  for 
the  approach  flow  of  the  laboratory  model  was  very  similar  to  the 
mean  value  of  the  land  and  water  profiles. 

Incoming  or  approach  wind  information  was  estimated  from 
synoptic  data  provided  by  Fleet  Numerical  Weather  Central,  the 
National  Meteorological  Center,  the  Japan  Meteorological  Agency, 
the  Royal  Observatory  at  Hong  Kong,  and  the  Naval  Weather  Service 
Environmental  Detachment  at  Cubi  Point.  This  data  base  was 
supplemented  by  the  wind  distribution  information  contained  in 
tropical  cyclone  warnings  issued  by  the  Joint  Typhoon  Warning 
Center,  Guam. 
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Figure  9.  Distribution  of 
synoptic  observations  from 
ships.  In  addition,  21 
observations  from  Cubi  Point 
recording  station  were 
incorporated.  Incoming  wind 
directions  range  from  225° 
to  255°. 
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Figure  10.  Normalized  profiles 
of  the  variation  of  wind  speed 
with  height  in  the  lowest 
boundary  layer  of  tropical 
cyclones  (after  Bates,  1978). 
The  dashed  line  indicates  the 
relationship  used  for  the 
present  study.  The  dotted  line 
indicates  the  normalized 
a  pp  ro  a  ch  -  fl  ow  profile  for  the 
1  aboratory  model  . 
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In  order  to  simplify  the  analysis  and  presentation  of 
results,  estimated  incoming  winds  (generally  ranged  from  25  to 
35  kt)  were  normalized  to  a  value  of  30  kt.  The  resulting 
harbor  winds  were  modified  proportionally  as  appropriate.  For 
example,  if  the  incoming  wind  was  35  kt  and  reduced  to  30  kt,  a 
ship  in  the  harbor  reporting  20  kt  would  be  proportionally  reduced 
to  17.1  kt.  Figure  11a  shows  the  results  of  the  synoptic  analyses 
of  composited  data  for  incoming  approach  winds  from  225°-255°. 

Figure  lib  presents  a  laboratory  model  composite  wind  distribu¬ 
tion  based  upon  the  incoming  winds  from  225°,  240°  and  255°  and 
drawn  for  the  comparable  11a  region. 

The  isotach  maximum  to  the  northwest  of  Cubi  Point  and  the 
primary  extension  of  the  maximum  to  the  east  are  evident  in  both 
Figures  11a  and  b.  A  secondary  extension  is  also  evident  to  the 
north  in  Figure  11a  but  is  more  pronounced  in  lib.  In  addition, 
the  decrease  in  winds  toward  the  western  shore  is  detectable  in 
both  figures.  While  these  pattern  similarities  agree,  the  synoptic 
isotach  values  are  generally  less  in  magnitude,  particularly  to 
the  north  of  the  maximum  and  in  the  sheltered  lee  locations  in  the 
eastern  section  of  the  harbor.  While  a  real  difference  might 
exist,  some  of  the  differences  might  be  attributable  to  over¬ 
estimates  of  the  wind  distribution  in  the  tropical  cyclone  warnings. 
An  overestimate  of  the  wind  distribution  would  contribute  to  an 
overestimate  of  the  incoming  wind.  Since  incoming  wind  and  harbor 
winds  were  proportionally  adjusted  to  an  incoming  wind  of  30  kt, 
an  overestimate  of  incoming  wind  would  lead  to  a  less  than  real 
harbor  wind.  In  an  examination  of  the  effects  of  30  typhoons 
crossing  the  Philippines,  Brand  and  Blelloch  (1973)  found  an 
average  surface  wind  decrease  of  33%.  If  the  effects  of  the 
Philippines  on  the  wind  and  wind  distribution  is  not  fully  taken 
into  account,  the  tropical  cyclone  warnings  containing  such 
information  might  be  overestimated. 

The  comparison  differences  in  the  sheltered  locations,  if  a 
real  difference  exists,  might  also  reflect  the  necessity  for  a 
higher  resolution  laboratory  model  using  output  from  the  larger 
scale  model  as  input. 
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igure  11.  (a)  Synoptic  analysis  of  composited  wind  data 

(141  observations)  normalized  to  30  kt  approach  wind  and 
70  ft  elevation;  analysis  is  only  over  region  where  data 
were  available,  and  is  for  the  approach  wind  range 
225°-255°.  (b)  Laboratory  model  composite  wind  distribu¬ 

tion  based  upon  the  incoming  winds  from  225°,  240°,  and 
255°;  analysis  is  only  over  the  region  shown  in  (a). 


4.  DYNAMIC  MODEL  OUTPUT 


A  diagnostic  examination  of  the  low-level  surface  wind  field 

was  also  made  by  means  of  a  one-level,  primitive  equation  model 

which  computes  mesoscale  influences  of  orography,  friction  and 

heating  on  surface  winds,  given  large-scale  synoptic  data.  The 

grid  interval  in  the  Subic  Bay  application  was  0.19  n  mi  and  the 

area  covered  was  a  4.3  x  4.3  n  mi  region  covering  the  eastern 

3 

section  of  the  bay. 

An  example  of  one  tropical  cyclone  case  that  was  diagnosti¬ 
cally  examined  is  discussed.  The  resulting  winds  are  for  an 
elevation  at  approximately  70  ft.  The  planetary  boundary  layer 

of  the  model  was  assumed  to  be  thoroughly  mixed  by  mechanical 

_  ? 

turbulence  and  the  drag  coefficients  used  were  1.3  x  10  for 

_  3 

land  and  1.1  x  10  for  water.  The  input  winds  at  the  top  of  the 
planetary  boundary  layer  were  250°  at  49  kt  with  resulting  low- 
level  winds  of  approximately  225°  at  30  kt.  Upper  air  tempera¬ 
tures  were  assumed  to  be  similar  to  the  mean  typhoon  temperatures 
of  Bell  and  Kar-siung  (1973). 

For  the  experimental  case  examined,  the  results  qualitatively 
compare  well  with  the  synoptic  and  laboratory  wind  distributions 
for  similar  low-level  approach  winds.  The  dynamic  model  results 
are  limited  to  the  eastern  sector  of  the  harbor  and  appear  to  show 
the  isotach  maximums  extending  toward  the  east  over  the  water 
regions  of  the  harbor  (Figure  12a).  The  comparable  synoptic 
composite  of  available  data  (Figure  12b)  also  shows  this  maximum 
to  the  east  north  of  Cubi  Point  runway  as  does  the  laboratory 
results  (Figure  12c).  The  wind  maximum  to  the  southwest  of  Cubi 
Point  runway  appears  also  to  extend  eastward  for  both  the  dynamic 
and  laboratory  distributions.  For  this  particular  wind  comparison, 
synoptic  data  was  not  available  to  the  southwest  of  Cubi  Point. 


3 

This  model  has  previously  been  applied  to  the  Juan  de  Fuca  and 
Georgia  Straits  in  British  Columbia,  using  a  5.4  n  mi  grid  size; 
see  Danard  (1977)  for  details. 
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Since  the  topography  to  the  west  of  the  bay  has  not  been 
considered  in  these  cases,  it  would  be  beneficial  to  run  a  larger 
scale  model  covering  the  whole  domain  and  then  use  the  output 
from  this  as  the  input  for  the  higher  resolution  grid. 
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5.  SUMMARY  AND  RECOMMENDATIONS 


An  attempt  has  been  made  to  examine  the  low-level  flow  of  a 
tropical  cyclone  in  the  vicinity  of  a  harbor  that  is  surrounded 
by  complex  terrain.  Several  approaches  were  used  to  deduce  the 
streamline-isotach  distributions  as  functions  of  incoming  wind 
directions.  Results  of  both  the  laboratory  experiment  and  the 
dynamic  model  show  promise.  However,  it  is  difficult  to  draw 
quantitative  conclusions  without  a  greater  number  of  more  detailed 
observations  from  both  inside  and  outside  the  harbor.  With  a 
detailed  observational  network,  the  Subic  Bay  area  would  provide 
a  setting  for  the  necessary  ground  truth,  especially  since  this 
major  Department  of  Defense  installation  is  threatened  by  an 
average  of  four  tropical  cyclones  each  year. 

The  mesoscale  (in  fact,  almost  microscale)  nature  of  the 
problem  is  certainly  challenging  from  the  viewpoint  that  very 
little  has  been  successfully  achieved  in  the  past  on  this 
diagnostic  scale,  although  this  is  precisely  the  scale  that 
seriously  concerns  many  local  forecasters,  particularly  under 
strong  wind  stress  conditions.  The  approaches  presented  in  this 
study  should  be  further  explored,  insofar  as  they  demonstrate 
potential  applications  to  the  many  other  harbors  worldwide  that 
are  affected  by  tropical  cyclones. 
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APPENDIX  A 

STREAMLINE-ISOTACH  DISTRIBUTIONS 

Streaml i ne-i sotach  distributions  (in  kt)  from  the  wind 
tunnel  laboratory  experiment  as  a  function  of  varying  wind  direc¬ 
tion  are  shown  in  Figures  A-l  through  A-7.  Wind  speed  values  are 
for  a  70  ft  elevation;  incoming  wind  speed  is  30  kt. 
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APPENDIX  B 

f 

STREAMLINE  AND  MAXIMUM  GUST  DISTRIBUTIONS  ' 


Streamline  and  maximum  gust  distributions  (in  kt)  from  the 
wind  tunnel  laboratory  experiment  as  a  function  of  varying  wind 
direction  are  shown  in  Figures  B-l  through  B-7.  Gust  values  are 
for  a  70  ft  elevation;  incoming  wind  speed  is  30  kt. 


